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Abstract 

This  effort  has  focused  on  verification  and  refinement  of  the  mechanism  transitions  at 
intermediate  temperatures  in  the  disk  alloy  Rene  104,  with  the  observation  of  microtwinning, 
continuous  faulting  and  dislocation  by-pass  at  successively  higher  temperatures  above  650°C. 
Evidence  for  the  twin  initiation  process  has  also  been  obtained  via  TEM  studies  of  specimens 
interrupted  after  small  strain  levels.  A  preliminary  model  for  the  novel  microtwinning  regime 
has  been  developed  that  appears  to  provide  reasonable  agreement  with  the  present  experimental 
creep  response  for  Rene  104  and  Rene  88.  Modeling  at  the  ab  initio,  atomistic  and  phase  field 
levels  is  providing  important  insight  into  the  activation  parameters  associated  with  the  observed 
deformation  mechanisms,  augmented  by  insight  from  3D  atom  probe  measurements  on  y  and  y' 
composition  and  ordering.  Single  crystals  of  Rene  104  have  successfully  been  grown,  and 
micro-tensile/compression  as  a  function  of  crystal  orientation  will  commence  in  the  next  year.  A 
novel  phase  field  model  of  directional  coarsening  (rafting)  during  high  temperature,  low  stress 
creep  of  blade  alloys  has  ben  developed.  This  model  accounts  for  the  local  stress  fields 
associated  with  matrix  dislocations  as  well  as  the  lattice  misfit,  and  demonstrates  promising 
qualitative  agreement  with  experiment.  Finally,  deformation  mechanisms  in  high  temperature 
creep  of  several  "generation  4"  single  crystal  blade  alloys  under  “rafting”  conditions  have  been 
studied  and  a  model  for  creep  under  these  conditions  developed. 

1.  Research  Objectives 

This  effort  comprises  a  coordinated  team  of  researchers  from  the  Ohio  State  University  (under 
Grant  #FA9550-05-l-0135),  University  of  Michigan  (under  Grant  #FA9550-05-l-0100)  and 
Johns  Hopkins  University  (under  Grant  #FA9550-05-l-0102).  The  goal  of  the  program  is  to 
develop  improved  models  that  will  (a)  incorporate  more  realistic  representation  of  the  relevant 
microstructures  and  micromechanisms,  (b)  enable  modeling  for  a  range  of  relevant  service 
conditions  (c)  address  time-dependent  deformation  in  both  disk  and  blade  alloys,  and  (d)  provide 
this  information  ev  to  the  component  design  process,  building  upon  the  paths  to  the  design 
process  created  in  the  DARPA  Accelerated  Insertion  of  Materials  (AIM)  program.  This  program 
is  coordinated  with  Pro  Tresa  Pollock  (University  of  Michigan)  is  leading  efforts  to  characterize 
blade  alloy  microstructures  and  is  providing  single  crystals  materials  for  the  rest  of  the  program 
and  Pro.  Kevin  Hemker  (Johns  Hopkins  Unviersity)  who  is  focusing  on  microtensile  tests  on 
single  crystal  variants  of  the  disk  alloy. 


2.  Accomplishments 

•  Completed  assessment  of  mechanisms  operative  during  creep  in  the  two-phase  superalloy  R104 
in  the  temperature  range  from  600-850°C.  As  a  function  of  temperature,  there  is  a  transition 
from  dislocation  activity  to  microtwinning,  extended/continuous  stacking  faults  and  dislocation 
climb  by-pass  at  the  highest  temperatures. 

•  Completed  detailed  analysis  of  microtwinning  mechanism,  which  is  the  dominant  mechanism 
at  the  upper  bound  temperature  for  service  of  disk  alloys  in  present  applications.  Microtwining, 
as  well  as  the  other  precipitate  shearing  modes,  is  controlled  by  diffusion-mediated  reordering. 

•  Determined  the  atomic-scale  processes  and  activation  barriers  associated  with  reordering  using 
the  nudged  elastic  band  method  coupled  with  first  principles  calculations  (VASP),  which 
indicate  that  the  diffusion  coefficient  for  reordering  should  be  similar  to  that  for  Ni-vacancy 
diffusion  in  the  ordered  Ll2  structure. 

•  Developed  a  phase  field  dislocation  dynamics  simulation  which  explains  the  transition  from 
dislocation  activity  to  microtwinning  as  a  function  of  precipitate  structure,  applied  stress  and 
fault  energies. 

•  Used  Z-contrast  imaging  to  detect  or  the  first  time  segregation  of  heavy  elements  to  microtwin 
and  stacking  faults  in  the  precipitates  after  creep  deformation.  The  fundamental  role  of  these 
heavy  elements  (C,  Mo  and  W)  in  modifying  these  fault  energies  has  been  explored  using  VASP 
calculations. 

•  Developed  a  model  for  high  temperature  creep  under  “rafting”  conditions  that  is  applicable  to 
“Generation  4”  type  single  crystal  superalloys. 

3.  Program  Summary  and  Present  Status 

During  this  period  of  the  program,  we  have  focused  on  improving  our  understanding  of  the 
mechanisms  of  creep  in  the  poly  crystalline  disk  superalloy,  Rene  104,  and  more  generally  on 
developing  the  modeling  tools  necessary  for  treating  critical  issues  in  the  high  temperature 
response  of  both  disk  and  blade  alloys.  These  activities  have  benefited  greatly  from  our 
interactions  with  industrial  colleagues  from  GE  Aircraft  and  Pratt- Whitney,  and  several 
collaborators  from  Wright  Laboratories. 

3.1  Mechanisms  of  Creep  Deformation  in  the  Disk  Alloy  Rene  104 

Creep  deformation  studies  have  been  conducted  on  the  advanced  disk  superalloy  Rene  104. 
Samples  for  creep  experimentation  were  extracted  from  a  fully  heat-treated  turbine  disk  supplied 
by  GEAE  and  P&W.  Previous  results  have  shown  that  microtwinning  is  the  principal 
deformation  mechanism  at  elevated  temperature  (>650°C),  which  corresponds  to  the  upper 
bound  of  service  temperatures  for  this  alloy  [1,2,5].  As  exemplified  in  Fig.  1,  twinning  proceeds 
by  the  passage  of  two,  closely-coupled  partial  dislocations  of  the  type  1/6<112>.  When  the 
partials  encounter  the  y'  precipitates,  high  energy  "complex  stacking  faults"  (CSFs)  are  created  in 
their  wake  due  to  the  formation  of  unfavorable  nearest- neighbor  bonds.  It  is  hypothesized  that 
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movement  of  the  twinning  partials  is  controlled  by  elimination  of  the  unfavorable  bonds  by  local 
atomic  rearrangements  (reordering  as  illustrated  in  Fig.  1)  and  is  responsible  for  the  temperature 
and  rate  dependence  exhibited  by  these  alloys.  Other  distinct  deformation  mechanisms  also 
reported  at  higher  temperature  involving  precipitate  shearing  may  also  be  controlled  by 
reordering  as  they  all  involve  shearing  of  the  y'  precipitates.  The  kinetics  associated  with 
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Fig  1:  (a)  HAADF  /STEM  image  of  two-layer  nanotwin  with  evidence  of  heavy 
element  segregation  (bright  columns  indicted  by  arrows)  and  (b)  hypothetical 
reordering-controlled  twinning  mechanism  in  the  y'  phase. 


reordering  have  been  explored  using  VASP,  and  initial  results  indicate  that  Ni-vacancy  diffusion 
and  reordering  should  have  similar  diffusion  coefficients  [13]. 


Several  heavy  elements  such  as  Mo,  Ta  and/or  W  are  added  to  commercial  alloys  since  they  tend 
to  improve  high  temperature  strength,  but  the  reason  for  this  improvement  has  not  been 
understood  in  detail.  Evidence  of  segregation  to  fault  and  microtwin  boundaries  has  also  been 
obtained  using  HAADF-STEM  imaging,  as  indicted  in  Fig.  1.  Interestingly,  segregation  is  only 
observed  in  the  y'  particles  [13].  Thus,  the  elements  may  be  swept  to  the  interfaces  in  the  course 
of  the  reordering,  ultimately  impeding  microtwin  thickening,  and  thereby  providing  enhanced 
creep  strength. 


3.2  Parametric  Studies  of  Mechanism  Transitions  Using  Phase  Field  Dislocation  Dynamics 


Insight  into  the  microstructural  factors  affecting  the  transition  from  matrix  dislocation  activity  to 
microtwinning  have  been  obtained  in  recent  experiments  coupled  with  phase  field  dislocation 
dynamics  modeling  [11],  We  have  observed  that  under  creep  (677°C  and  724MPa  and 
microstructure  (fast  cooling/fine  microstructure)  conditions  for  which  microtwining  is  the 
dominant  process  at  2%  strain,  a  distinctly  different  substructure  is  observed  at  smaller  strain 
levels.  Detailed  dislocation  indicates  that  1/2<110>  matrix  dislocations,  all  having  the  same 
Burgers  vector  in  a  given  grain,  dissociate  widely  into  Shockley  partial  dislocations  that  have  a 
Burgers  vector  of  the  type  a/6<l  12>.  The  leading  partial  dislocations  appear  to  have  "percolated" 
between  the  secondary  y'  precipitates,  leaving  them  unsheared.  While  the  process  whereby  this 
structure  subsequently  evolves  into  microtwins  is  presently  not  clear,  the  decorrelation  of  Vi [110] 
matrix  dislocations  nevertheless  appears  to  be  a  necessary  precursor  to  microtwinning. 

The  interaction  of  a/2<110>  dislocation  with  precipitates  and  the  various  kinetic  pathways 
associated  with  their  dissociation  have  been  explored  using  the  so-called  microscopic  phase-field 
(MPF)  model.  The  MPF  enables  treatment  of  complex  y'  shearing  processes  without  any  a  priori 
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assumptions  about  dislocation  geometry,  core  structure,  and  formation  of  stacking  faults.  The 
result  of  one  such  simulation  for  a  single  screw-oriented  dislocation  as  it  interacts  with  a 
secondary  y'  precipitate  population  with  a  size  (diameter)  of  300  nm  and  an  interparticle  spacing 
of  75nm  is  shown  in  Figure  2. 

As  the  applied  stress  in  this  case  is  below  the  Orowan  stress  for  the  perfect  dislocation  (AB),  it  is 
sufficient  to  drive  the  leading  partial  (A6)  into  the  channel  because  of  the  smaller  line  tension 
(which  is  approximately  proportional  to  b2’,  where  b  is  the  Burgers  vector  of  the  partial).  Note  the 
stress  resolved  on  the  trailing  partial  is  only  half  that  on  the  leading  partial  for  the  orientation 
chosen.  As  shown  in  Fig.  2a,  the  leading  Shockley  sweeps  out  and  creates  an  intrinsic  stacking 
fault  in  the  matrix  while  the  trailing  partial  remains  at  the  entrance  to  the  y  channel.  For  the 
same  microstructure,  but  increasing  the  matrix  stacking  fault  energy  10  fold  to  100mJ/m2,  leads 
to  a  quite  different  result,  as  shown  in  Fig.  2b.  Under  these  circumstances,  when  the  dislocation 
makes  its  way  to  the  confines  of  the  y  channel,  due  to  the  increase  in  matrix  stacking  fault 
energy,  the  dislocation  does  not  dissociate  into  Shockley  partials.  Instead,  the  dislocation 
remains  intact  as  a  full  a/2<l  10>  and  is  immobilized. 

These  preliminary  results  of  the  phase  field  model  demonstrate  several  important  factors  that 
trigger  the  decorrelation  of  the  Vi  [110]  dislocations,  and  leading  ultimately  to  deformation  via 
microtwining.  These  factors  include  the  matrix  stacking  fault  energy,  as  high  fault  energy 
inhibits  decorelation,  and  ultimately  microtwinning.  At  a  stress  below  the  Orowan  bowing  stress 
for  Vi  [110]  dislocations,  decorrelation  will  be  favored  for  microstructures  in  which  the  spacing 
between  secondary  particles  is  smaller.  If  the  stress  is  in  excess  of  the  Orowan  bowing  stress, 
then  the  matrix  dislocations  remains  undissociated,  regardless  of  stacking  fault  energy.  Note 
also  that  the  absence  of  twins  at  this  higher  stress  level  would  be  difficult  to  understand  if 
dislocation  activity  and  microtwinning  were  operative  in  parallel.  The  phase  field  simulations 
instead  suggest  that  if  the  microstructure  and  stress  conditions  are  favorable  to  Vi  [110]  matrix 
activity,  rather  than  decorrelation,  then  microtwining  will  not  be  operative.  The  orientation  of 
the  applied  shear  stress  relative  to  the  operative  Burgers  vector  should  have  a  significant  role  on 
decorrelation.  This  indicates  that  decorrelation  and  microtwinning  may  only  be  favored  over  an 
relatively  narrow  range  of  applied  shear  stress  orientation.  These  attributes  will  be  incorporated 
into  unique  constitutive  formulations  that  are  being  developed  for  the  crystal  plasticity  modeling 
by  Ghosh. 

3,3  Deformation  Mechanisms  and  Microstructure  Evolution  for  Rafted  Microstructures 

Our  modeling  efforts  have  been  coordinated  with  deformation  studies  at  University  of  Michigan 
by  Prof.  Tresa  Pollock,  also  supported  by  the  MEANS2  program.  As  is  typical  of  the  single 
crystal  blade  alloys,  directional  coarsening  of  the  y’  precipitates  occurs  during  creep  exposure. 
Dislocation  substructures  are  being  studied  at  various  stages  of  creep  deformation  to  obtain 
critical  information  on  rate-limiting  mechanisms  in  the  rafting  regime.  A  ground-breaking  phase 
field  model  that  reproduces  these  remarkable  morphological  changes  has  also  been  developed 
under  this  program  in  collaboration  with  Prof.  Yunzhi  Wang  and  his  graduate  student,  Ning 
Zhou,  and  initial  results  of  this  model  are  described  below  and  in  [6,8].  The  simulation 
technique  is  based  on  an  integrated  phase  field  model  that  for  the  first  time  incorporates 
simultaneously  spatiotemporal  evolution  of  both  precipitate  morphology  and  dislocation 
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structures.  The  initial  configurations  consisting  of  cuboidal  y'  particles  and  dislocations  in  y 
channels  are  constructed  according  to  experimental  observations  and  phase  field  simulation  of 
dislocation  filling  process  in  the  y  channels.  The  spatial  variation  of  chemical  potential  of  solute 
atoms  is  evaluated  based  on  local  concentration  and  stress  and  diffusion  fluxes  in  different  y 
channels  are  analyzed. 


ISF  Energy 
100  mJ/m2 


Figure  2.  Phase  field  simulation  of 
a/2<110>  dislocation  (AB)  interaction 
with  a  secondary  y'  precipitate  at  an 
applied  shear  stress  of  115MPa,  300nm  y' 
precipitate  size,  and  75nm  y  channel  width 
spacing.  The  simulation  shows  the  effect 
of  tertiary  y'  precipitates  and  stacking  fault 
energy  a)  10mJ/m2  ISF  and  b)  100mJ/m2 
ISF  energy  w/tertiary  y'.  Periodic 
boundary  condition  is  applied  in  the 
vertical  direction.  Note  the  spontaneous 
dissociation  under  stress  of  the  matrix 
dislocation  as  it  squeezes  through  a  narrow 
matrix  channel. 


In  this  model,  dislocation  motion  in  y  channels  is  driven  by  the  local  stress.  This  local  stress 
includes  contributions  from  lattice  misfit,  the  external  load,  and  misfit  dislocations  at  the  y/y' 
interface.  Simultaneously,  the  diffusion-mediated  evolution  of  microstructure  is  determined  by 
both  chemical  driving  forces  and  elastic  interactions  with  dislocations,  as  well  as  the  effect  of 
lattice  misfit.  Through  this  integrated  approach,  rafting  of  the  y1  phase  has  been  predicted  under 
various  conditions  of  misfit  sign  and  loading  type,  and  the  results  agree  with  our  experiments  and 
other  analyses.  Some  results  are  shown  in  Figs.  3  and  4.  It  has  also  been  found  that  modulus 
mismatch  is  not  a  necessary  condition  for  rafting,  as  has  been  assumed  previously  in  the 
elasticity-based  models  for  rafting.  Under  the  homogeneous  modulus  assumption  (ignoring  the 
modulus  mismatch  between  y  and  y'  phase),  the  simulation  results  show  rafting  of  y'  precipitates 
is  closely  related  to  the  type  and  distribution  of  dislocations  at  the  y/y'  interfaces.  The  sign  of 
lattice  misfit  and  the  direction  of  the  applied  load  determine  the  dislocation  configurations  in  the 
y-channels  and  the  rafting  directions.  It  is  found  that  dislocations  prefer  normal  channels  (normal 
to  the  direction  of  applied  load)  for  a  positive  misfit  alloy  under  compression,  and  a  negative 
misfit  alloy  under  tension  which  later  develop  N  type  rafting,  while  parallel  channels  were 
preferred  by  dislocations  in  a  positive  misfit  alloy  under  tension  and  a  negative  misfit  alloy  under 
compression  which  develop  P  type  rafting. 
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Figure  3:  Phase  field  model  result  showing  N  type  rafting  for  alloys  with  -0.3%  lattice  misfit 
under  152  MPa  tensile  stress  along  [001]  after  11  hours. 


Figure  4:  Simulation  for  P  type  rafting  for  alloys  with  +0.3%  lattice  misfit  under  152  MPa 
tensile  stress  along  [001]  after  7  hours. 


The  model  has  further  been  developed  to  include  modulus  mismatch,  using  an  effective 
transformation  strain  to  describe  the  channel  plasticity.  A  parametric  phase  field  simulation 
study  [12,14]  has  resulted  in  the  following  conclusions: 

•  When  both  channel  plasticity  and  modulus  mismatch  are  considered,  the  type  of  rafting  is 
determined  solely  by  the  direction  of  the  applied  stress  and  the  sign  of  the  lattice  misfit  (i.e.,  the 

sign  of  — ).  The  sign  of  the  modulus  mismatch  does  not  change  the  rafting  type,  which  is  true 
<5 

even  when  the  difference  in  the  shear  modulus  is  about  15%,  i.e.,  (C"-C,J  ,85%- 

•  The  P  type  rafting  completes  faster  than  the  N  type  rafting.  The  reason  is  that  for  the  P  type 
rafting,  solute  atoms  flow  from  two  vertical  channels  to  one  horizontal  channel  leading  to  a  rod; 
while  for  the  N  type  rafting,  solute  atoms  flow  from  one  horizontal  channel  to  two  vertical 
channels  leading  to  a  plate.  In  other  words,  forming  a  plate  takes  longer  time  than  forming  a  rod. 

•  Modulus  mismatch  between  precipitate  and  matrix  does  alter  the  rafting  kinetics.  If  the 
precipitate  phase  has  higher  elastic  stiffness  than  that  of  the  matrix,  the  rafting  process  will  take 
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less  time  to  complete  as  compared  to  a  system  having  softer  precipitates  (assuming  the  same 
applied  stress,  misfit  and  dislocation  density). 

•  Dislocation  density  will  also  affect  the  kinetics  of  rafting.  The  higher  the  channel  dislocation 
density,  the  faster  the  rafting  process.  The  density  of  dislocations  is  determined  by  the  magnitude 
of  the  applied  stress,  lattice  misfit  and  the  starting  microstructure  (y'  particle  size  and  y  channel 
width). 

3.4  Creep  Deformation  of  Advanced  Single  Crystals  in  the  Rafting  Regime 

The  alloys  studied  cover  a  range  of  composition  representative  of  advanced  “Generation  4” 
single  crystals.  In  these  alloys  directional  coarsening  of  the  y’  precipitates  occurs  during  creep 
experiments  conducted  at  950°C  and  295  MPa.  Dislocation  substructures  are  being  studied  at 

various  stages  of  creep 
deformation  in  sections  parallel  to 
the  applied  tensile  stress  to  obtain 
critical  information  on  rate- 
limiting  mechanisms  in  the  rafting 
regime.  No  comprehensive  creep 
models  have  been  developed  for 
this  stress  /  temperature  regime, 
which  has  received  relatively  limited  attention  beyond  studies  that  address  the  directional 
coarsening  phenomenon.  With  regard  to  the  rate-limiting  mechanisms  three  possibilities  are 
being  considered:  (1)  creep  controlled  by  viscous  glide  and  climb  of  dislocations  across  the  y’ 
rafts;  (2)  creep  controlled  by  viscous  glide  and  climb  in  the  matrix  y  phase  and  (3)  creep 
controlled  by  diffusion-limited  reactions  at  y  -  y’  interface  that  produce  dislocations  for  shearing 
of  the  y’  rafts. 


Alloy 

Ni 

Al 

Ta 

Re 

W 

Co 

Cr 

UM 

F20 

67 

13.9 

3.5 

2.8 

1.5 

1.0 

2.5 

8.0 

UM 

F30 

59 

13.8 

3.5 

2.7 

1.5 

1.0 

10.5 

8.0 

Table  1:  Compositions  for  two  experimental  "generation 
4"  alloys  under  investigation. 


6x10" 

5x1 

Q  3x10" 

J  2x10" 

<u 

8  1x10'* 

cn 

Q  n 


i  ’  i 


i  ■  i  ■  > 


Gliding  dislication  in  y  matrix 
In  rafted  y'  phase 


UMF30 


I 


UMF20 

z 


The  most  recent  studies  have  focused  on  two  alloys 
designated  UM-F20  and  UM  F-30,  whose 
compositions  are  similar  except  for  the  Co  level  and 
are  listed  in  Table  1.  We  note  the  presence  of 
dislocations  in  the  y’  rafts,  at  the  y  -  y’  interfaces  and 
in  the  gamma  matrix.  The  density  of  dislocations 
that  appear  to  be  mobile  in  each  of  the  phases  was 
measured,  and  the  results  are  shown  in  Figure  5. 
Note  that  with  changes  in  minimum  creep  rate,  there 
are  changes  in  the  matrix  dislocation  density,  but  no 
changes  in  the  density  within  the  precipitates. 
Furthermore,  the  creep  rate  decreases  in  UMF30 
compared  to  UMF20,  with  increasing  dislocation 
density  in  the  matrix  and  also  with  increasing  Co  in 
the  matrix.  This  seems  to  indicate  that  dislocations 
accumulate  in  the  matrix  until  favorable  conditions  develop  at  the  interface  for  pairing  of 
dislocations  in  a  configuration  that  will  result  in  shearing  of  the  y’  rafts.  Based  on  the 
experimental  observations  reported  here,  two  possibilities  for  the  rate  controlling  mechanism 
during  high  temperature  creep  of  superalloys  in  the  rafting  regime  can  be  envisaged:  (1)  the 


0  1x10*  2x10*  3x10*  4x10*  5x10*  6x10* 

Experimental  Minimum  Creep  Rate  (Is) 

Figure  5:  Measured  dislocation 
densities  in  the  y  and  y'  phases  after 
creep  in  rafted  microstructures. 
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viscous  motion  of  dislocation  pairs  through  the  y’  rafts  and  (2)  periodic  (time  dependent) 
penetration  of  dislocations  from  the  interfacial  networks  into  the  rafted  precipitates.  Srinivasan 
and  coworkers  have  suggested  that  pairwise  cutting  can  occur  within  ’  rafts  by  a  combined  climb 
and  glide  mechanism  in  CMSX-6  at  1020°C  and  80  MPa  [44],  The  creep  rate  is  then  limited  by 
the  viscous  motion  of  the  climbing  dislocation  through  the  raft.  In  this  case,  the  experimental 
creep  rate  coincided  well  with  the  calculated  strain  rate,  based  on  a  model  of  the  viscous 
shearing.  A  similar  calculation  was  performed  for  the  present  Ru-containing  experimental 
superalloys  using  the  climb  velocity  determined  by  Srinivasan  and  coworkers: 

y  =  0.5pbvc 

where  is  the  dislocation  density,  b  is  the  Burgers  vector,  and  vc  is  the  climb  velocity  determined 
for  CMSX-6  as  4x1  O'8  m/s.  For  the  present  Ru-containing  alloys,  the  climb  velocity  is  assumed  to 
be  similar  based  on  a  lower  creep  temperature  but  higher  magnitude  of  the  externally  applied 
stress.  The  density  of  dislocations  within  the  rafts,  ,  was  determined  from  TEM  images  taken 
longitudinally  in  the  gage  section  of  UM-F20  after  creep  rupture,  and  was  determined  to  be 
3xl012  /m2.  Therefore,  the  creep  rate  is  estimated  as  1.510'5  /s.  The  experimentally  determined 
creep  rate  for  UM-F20,  however,  was  several  orders  of  magnitude  lower,  at  6xl0'9  /s. 

The  considerable  difference  between  the  calculated  and  experimentally-determined  creep  rate 
suggests  that  the  viscous  motion  of  dislocation  pairs  through  the  y’  rafts  is  not  the  rate  limiting 
process  for  this  alloy  at  950°C  and  290  MPa.  The  other  possibility  for  the  rate  limiting  creep 
mechanism  is  the  time  dependent  penetration  of  the  dislocations  through  the  interfacial 
dislocation  networks  and  into  the  rafts.  The  time  dependence  may  be  a  result  of  either  (a)  the 
frequency  at  which  two  parallel  a/2<110>  dislocations  arrive  at  the  y’  interface,  or  (b)  the 
frequency  at  which  a  gliding  matrix  a/2<110>  dislocation  encounters  a  compatible  a/2<110> 
network  dislocation.  The  former  is  unlikely  based  on  the  present  TEM  observations  that  show 
excess  dislocations  present  in  the  matrix  channels.  Rather,  the  present  microscopy  studies 
suggest  that  the  time  dependence  of  penetrating  pairs  of  dislocations  is  the  most  likely 
controlling  mechanism  for  creep  in  this  class  of  alloys  at  the  conditions  evaluated.  The  network 
structure  is  rapidly  established  along  the  y’  interfaces  during  primary  creep;  further  entry  of 
dislocations  that  sustain  the  creep  process  would  elicit  interactions  between  the  existing  network 
structure  and  the  gliding  matrix  dislocations.  Considering  the  deposition  of  a  dislocation 
segment  by  glide  of  a  matrix  a/2<110>,  it  is  likely  that  its  position  at  the  interface  will  lie 
between  two  pre-existing  network  dislocations.  If  the  dislocation  is  of  the  correct  (similar)  sign, 
climb  of  this  deposited  dislocation  onto  an  adjacent  glide  plane  with  a  network  segment  would 
provide  the  necessary  pairwise-cutting  source  for  raft  penetration.  Note  that  significant  climb 
has  already  occurred  within  the  network  to  achieve  a  desirable  misfit-reducing  configuration  of 
the  net.  Another  local  climb  event  of  significance  would  involve  dislocations  of  opposite  signs 
climbing  to  an  adjacent  plane  in  order  to  annihilate  the  network  dislocation,  resulting  in  a 
reduction  in  dislocation  density  and  localized  coarsening  of  the  network.  Because  dislocations 
entering  a  channel  with  an  established  network  have  an  equal  chance  of  being  of  similar  or 
opposite  sign  to  neighboring  dislocations,  the  rate  of  production  of  a  dislocation  pairs  that  can 
shear  rafts  would  be  directly  related  to  the  rate  of  coarsening  of  an  established  network  during 
creep. 


8 


Considering  that  the  creep  rate  is  limited  by  the  climb  velocity  of  dislocations  into  critical 
shearing  configurations  at  the  scale  of  the  interfacial  network,  then  the  overall  creep  rate  can  be 
estimated  from  the  following: 


£ 


where  is  the  interfacial  dislocation  density  and  dx/dt  represents  the  average  glide  distance  per 
raft  shearing  event  per  unit  time.  The  dx/dt  term  will  be  governed  by  the  climb  distance,  A  and 
climb  rate,  vc,  of  the  dislocations  at  the  y’  interface,  as  well  as  the  width  of  the  raft,  (3.  Thus: 

dx  p 

- Vc 

dt  A 


with  climb  velocity  as 

DQ.ab 

Vc  = - ^ - 

frkT 

where  D  is  the  diffusivity,  is  the  atomic  volume,  b3,  k  is  Boltzmann’s  constant,  and  T  is  the 
temperature  [45] .  Thus  the  creep  rate  would  then  be 

kT  A. 


It  is  necessary  to  consider  only  the  density  of  dislocations  involved  in  generating  shearing  pairs 
of  dislocations  from  the  interfacial  network.  Thus  only  a  fraction  of  the  dislocations  in  the 
vetcoopK  will  experience  an  interaction  with  the  appropriate  type  of  gliding  dislocation;  this 
fraction  is  designated  by  a  pre-exponential  factor,  Z, 


ZfipDotd 

m 


Recent  interdiffusion  experiments  on  Ni-Re  binary  alloys  suggest  that  the  interdiffusion 
coefficient  at  950°C  is  approximately  D  =  6xl0'18  m2/s  [46].  The  density  of  interfacial 
dislocations  is  approximately  1015/m2,  while  (3/A  =  10  for  the  present  alloys.  The  stress  is 
approximately  290  MPa  in  the  present  experiments.  Considering  a  creep  rate  of  approximately 
10"8  /s,  an  estimate  of  the  pre-exponential  factor,  Z,  on  the  order  of  10"6  is  obtained.  Not 
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surprisingly,  this  suggests  that  only  a  few  segments  of  the  dislocations  within  the  interfacial 
network  are  involved  in  the  accumulation  of  strain  during  high  temperature  creep.  So  it  is  worth 
considering  the  spacing  between  such  active  links  in  the  network  present  with  density  act.  With 
Z  =  10'6,  approximately  one  dislocation  out  of  every  106  of  the  network  dislocations  is  active. 
This  spacing  between  active  links  should  be  of  the  same  order  of  magnitude  as  observed  in 
Figure  6  and  the  spacing  that  is  calculated  from  the  density  of  dislocation  pairs  in  the  y’  rafts, 
10 12  /m2,  corresponds  to  approximately  1  dislocation  pair  shearing  the  raft  per  1  //m.  The 
fraction  of  active  dislocations  was  determined  from  the  experimental  spacing  between  shearing 
events  and  the  interfacial  dislocation  spacing  as  approximately  5*10'4,  quite  similar  to  the 
calculated  value,  Z,  which  was  approximately  10'6.  Thus,  this  model  is  physically  reasonable. 
One  implication  is  that  a  high  resistance  of  the  y’  phase  to  shearing  is  important  to  improving 
creep  properties.  An  alloy  that  has  a  high  precipitate  shearing  resistance  will  accumulate  excess 
matrix  dislocations.  A  high  precipitate  shearing  resistance  will,  in  turn,  be  dependent  on 
precipitate  composition  and  the  associated  APB  and/or  complex  fault  energies.  Finally,  creep 
resistance  can  be  improved  with  the  addition  of  elements  with  slow  interdiffusion  coefficients  in 
the  matrix,  in  order  to  slow  the  climb  kinetics  in  the  near  interfacial  regime.  Rhenium  additions 
are  well  known  to  be  effective  in  this  regard. 
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